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Abstract 

Nonspherical perturbation theory has been necessary to understand the 
meaning of radiation in spacetimes generated through fuhy nonhnear nu- 
merical relativity. Recently, perturbation techniques have been found to be 
successful for the time evolution of initial data found by nonlinear methods. 
Anticipating that such an approach will prove useful in a variety of problems, 
we give here both the practical steps, and a discussion of the underlying the- 
ory, for taking numerically generated data on an initial hypersurface as initial 
value data and extracting data that can be considered to be nonspherical 
perturbations. 

I. INTRODUCTION 

The formation of a black hole is, in principle, one of the most efficient mechanisms for 
generation of gravitational waves. Such sources tie together two major research initiatives. 
Laser interferometric gravity wave detectors |^ hold out a promise of the detection of grav- 
itational waves from astrophysical events. To interpret the results of the gravitational wave 
signals, and to help find signals in the detector noise, a broad and detailed knowledge will 
be needed of astrophysical gravitational waveforms. This is one of the underlying motiva- 
tions for the "grand challenge" |^ in high performance computing, aimed at computing the 
coalescence of black hole binaries. 

Evolving numerical spacetimes and extracting outgoing radiation waveforms is indeed 
a challenge. In a straightforward numerical approach, a good estimate of the asymptotic 
waveform requires long numerical evolutions so that the emitted waves can be propagated 
far from the source. The necessary long evolutions are difficult for a number of reasons. 
General difficulties include throat stretching when black holes form, numerical instabilities 



associated with curvilinear coordinate systems, and the effects of outer boundary conditions 
which are approximate. 0] 

We suggest here that at least part of the cure for this problem may lie in the use of 
the theory and techniques of nonspherical perturbations of the Schwarzschild spacetime 
( "NPS"). By this we mean the techniques for treating spacetimes as deviations, first order in 
some smallness parameter, from the Schwarzschild spacetime. These techniques differ from 
"linearized theory" which treats perturbations of the spacetime from Minkowski spacetime 
and which cannot describe black holes. The basic ideas and methods were set down by many 
authors and lead to "wave equations" for the even parity |^ and odd parity perturbations. 

NPS has been used to compute outgoing radiation waveforms from a wide variety of black 
hole processes, including the scattering of waves P] , particles falling into a hole , and stellar 
collapse to form a hole P). The general scheme of NPS also underlies the techniques for 
extraction of radiation from numerically evolved spacetimes 0. NPS computations have 
recently been used in conjunction with fully numerical evolution, as a code test |TD[ and as 
a strong- field radiation extraction procedure 0. 

Here we are interested in another sort of application of NPS theory. To understand such 
applications we consider an example: Two very relativistic neutron stars falling into each 
other, coalescing and forming a horizon, as depicted in Fig. 1. The curve "hyper surf ace," 
in Fig. 1, indicates a spacelike "initial" surface. The spacetime can be divided into three 
regions by this initial surface and the horizon. The early evolution, in region I, below the 
initial hypersurface, is highly dynamical and nonspherical. Spherical perturbation theory is 
clearly inapplicable. Above the initial surface the spacetime remains highly nonspherical in 
region II inside the event horizon, but outside the event horizon, in region III, it may be 
justified to consider the spacetime to be a perturbation of a Schwarzschild spacetime. This 
is essentially guaranteed if the initial hypersurface is chosen late enough, in some sense, after 
the formation of the horizon. The evolution in region III, then, is determined by cauchy 
data on the initial hypersurface exterior to the horizon. It is important to note that this 
is made possible by the fact that the horizon is a causal boundary which shields the outer 
region from the dynamics of the highly nonspherical central region. 

The scheme inherent in this division of spacetime has the potential greatly to increase the 
efficiency of the computation of the radiation generated when strong field sources form black 
holes. If one starts from the cauchy data on the initial hypersurface, one can evolve forward 
in time with the linear equations of perturbation theory. Many of the long-time evolution 
problems of numerical relativity are avoided and the interpretation of the computed fields 
in terms of radiation is immediate. 

The approach suggested would then seem to be: Use numerical relativity up to the initial 
hypersurface; use the techniques of nonspherical perturbations in the future of the initial 
hypersurface. In fact, the efficiency that can be achieved may be even greater. In the early, 
highly nonspherical, pre- initial hypersurface phase of the development of the spacetime, there 
may be relatively little generation of gravitational radiation. By using a computational tech- 
nique which suppresses the radiative degrees of freedom one may be able to compute the 
early stages of evolution relatively easily. There are two very recent examples of just such 
applications of this viewpoint. Price and PuUin |1T[ used as initial data the Misner's |T^ 
solution to the initial value equations for two momentarily stationary black holes. Abra- 
hams and Cook [O considered two holes moving towards each other, and used numerical 
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values of the initial value equations. In neither case was there any use of fully nonlinear nu- 
merical evolution. The rather remarkable success of both computations suggests that there 
is something robust about the underlying idea of separating horizon-forming astrophysical 
scenarios into an early phase with no radiation and a late phase with small deviations from 
sphericity outside the horizon. It is plausible that the bulk of the radiation in most processes 
is generated only in the very strong-field interactions around the time of horizon formation 
and that radiation generation in the early dynamics can be ignored. One would, however, 
think that strong radiation would be emitted during the stages at which the early horizon 
is very nonspherical and at which time nonspherical perturbation theory would seem to be 
inapplicable. There should be a tendency for this "early" radiation, produced very close 
to the horizon, to go inward into the developing black hole, so that the application of non- 
spherical perturbation theory to the exterior really requires that on the initial spacetime 
the perturbation are small only well outside the horizon. It would seem that something 
of this sort would have to be happening to explain the accuracy of the Price-PuUin and 
Abrahams-Cook results. 

Whether or not many problems can be treated with no use of fully numerical evolution, it 
appears clear to us that these perturbation methods will be applied to a variety of problems 
in which data on the initial hypersurface is available numerically. The primary purpose 
of this paper is to provide justification and background for earlier work on this subject 
and a clear recipe for future applications. In the next section we discuss the meaning, and 
limitations, of extracting a "perturbation" from this numerical data and computing radiated 
energies. The explicit process of extracting the perturbations from the numerical data is 
given in Sec. III. In Sec. IV we demonstrate the use of this procedure via application to a 
specific example, the Misner initial data. 



II. INITIAL DATA AS SCHWARZSCHILD PERTURBATIONS 

We outline here the formalism for perturbation theory based on work by Regge and 
Wheeler ^ and by Zerilli 0] , but we will draw heavily on the gauge invariant reformulation 
of those earlier works by Moncrief WM. Our starting point is an initial hypersurface which 



can be taken as a surface of constant Schwarzschild time. We assume that the coordinates 
on that surface are almost Schwarzschild coordinates r, 6, and we assume that the 
values are known, on this hypersurface and in these coordinates, for the 3-metric jij and the 
extrinsic curvature Kij. The conditions for finding such a hypersurface and such coordinates 
will be made explicit in Sec. III. 

Underlying perturbation theory is the idea of a family of metric functions (7^,y(x";e), 
depending on the parameter e, which satisfy the Einstein equations for all e, and which, 
in the limit e ^ 0, become the Schwarzschild metric functions, such as Qrr = S^^. (Here 
5* = 1 — 2M/ r and M is the mass of the Schwarzschild spacetime; we use units throughout 
in which c = G = 1.) NPS theory amounts to the approximation 

d 

g^,u{x°'■, e) ~ gfiuix""] e)|,=o + e—g^^,{x"; e)|,=o . (1) 
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A. Choice of expansion parameter 



It is of some practical importance to realize that the choice of the expansion parameter 
can have a considerable effect on the range over which perturbation theory gives a good 
approximation. Let us imagine that we introduce a new parameter e' which is a function of 
e such that de' / de approaches unity as e — 0. If we take e' to be the basis of our perturbation 
approach, the approximation becomes 

g{x"- e) = g{x^- e(e')) = g,u{x^; e(e'))|.'=o + e'^g,,{x-- e(e'))|e'=o + 0{e'^) 

= g,,{x^- e(e'))|.=o + [t^9,.{x^; e(e'))|.=o] } + 0{e'^) . (2) 

At e = the derivative of g^^, with respect to e and with respect to e' have the same values, 
so for a given spacetime — that is, for a given value of e — the nonspherical perturba- 
tion in d^) differs from that in ([l|) by the factor {e'/e}. Computed energies (which are 
quadratic in the nonspherical perturbations) will differ by the square of this ratio. Different 
choices of parameterization will change this factor and affect the accuracy of the linearized 
approximation. 

To show the effects of this parameterization dependence, we take as an example Misner 
data inil for two holes. The initial separation of the holes, in units of the mass of the 



spacetime, is described by Misner's parameter /xq- The metric perturbations, however, are 
not analytic in /xq as /xq 0, so Hq cannot be used as the expansion parameter in (|l|). The 
actual expansion parameter used by Price and Pullin, was a function of /Xq denoted ^2- We 
consider here what would be the results of perturbation theory done with the expansion 
parameter 

(3) 



1 — kK2 



Figure 2 shows the results, along with the energies computed by numerical relativity ap- 
plied to full nonlinear evolution |15[ . For all choices of k the agreement between perturbation 



theory and numerical relativity is good at sufficiently small initial separation (sufficiently 
small /io); but as /io grows larger, the agreement increasingly depends on the which param- 
eterization is used. The = parameterization, the parameter of the Price-PuUin paper, 
is a reasonably good approximation even up to separations (/io > 1.36) for which the initial 
apparent horizon consists of two disjoint parts. For positive values of k the agreement is 
less impressive, while for k = —4, it appears that perturbation theory is giving excellent 
answers for initial data that are very nonspherical. Clearly the k = —4 parameterization 
is "better," at least for the purpose of computing radiated energy. There exist yet better 
choices; in principle a parameterization could be found for which the energy computed by 
linearized theory is perfect for any initial separation. The crucial point is that we have no 
a priori way of choosing what is and what is not a good parameterization. The choice of 
expansion parameter K2 was made in the Price-PuUin analysis, because it occurred naturally 
in the mathematical expressions for the initial geometry. There was no a priori reason for 
believing it to be a particularly good, or particularly bad parameterization. This point will 
be discussed again, in connection with numerical results presented in Sec. IV. 

The fact, demonstrated in Fig. 2, that the results of linear perturbation theory are arbi- 
trary may seem to suggest that perturbation answers, from a formal expansion or numerical 
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initial data, are of little value. It should be realized, however, that the arbitrariness exhib- 
ited in Fig. 2 is simply a demonstration of the fact that linearized perturbation results are 
uncertain to second order in the expansion parameter. The fact that the results for different 
parameterizations start to differ from each other around fiQ ~ 1.5 simply signals that K2 is 
around unity. (In fact, k,2 ~ 0.24 for fiQ = 1.2.) Higher order uncertainty is an unavoid- 
able feature in the range where the expansion parameter is of order unity. But there is a 
potential misunderstanding about the meaning of "expansion parameter around unity." To 
see this consider a change to a new expansion parameter e = 10~^ * K2. The new expansion 
parameter e is of order unity for hq ^ 7, yet we know that perturbation fails dramatically 
for such a large value of fiQ. The issue here is that we need some way of ascribing an appro- 
priate "normalization" to the expansion parameter. A sign that the normalization is good 
is that physically-based measures of distortion start getting large for e around unity. If we 
had reliable measures of this type then we could have some confidence about the range of 
the the expansion parameter for which we could neglect second order uncertainty, whether 
due to parameter arbitrariness or the omission of higher order terms in the calculation. 
One can formulate interesting measures for the normalization of the expansion parameter, 
such as the extent to which the linearized initial conditions violates the exact Hamiltonian 
constraint [^. Most such measures are useful only for finding a very rough normalization 
for K2 (equivalently, for roughly finding the range in which linearized perturbation theory 
is reliable). The only reliable procedure for this is to carry out computations of radiated 
waveforms and energy to second order in the expansion parameter. The ratio of second 
order corrections to first order results gives the only direct measure of the reliability of per- 
turbation results. If one computes an energy for which the second order correction to the 
first order result is 10%, then one knows that the third order correction (due to a change 
in parameterization or an inclusion of third order terms in the computation) will be on the 
order of 1%. 



B. Treating nonlinear initial data as a perturbation expansion 

We turn now to the central question of this paper: How does one apply perturbation 
theory to numerically generated initial data? To do this we consider our numerical initial 
data to be initial data for a solution in a parameterized family g^i,{x'^] e) corresponding to 
^ — Cnum- The application of perturbation theory is equivalent to replacing g^yix"^] Cnum )by 

d 

g^y (x" ; e) | e=o + -Q-^9^iu (a;" ; e) | ,=oenum • (4) 

An added familiar complication is that we can introduce a family of coordinate transfor- 
mations = x°(a;^ ; e) which reduces to = x'^ for e ^ 0. Such a transformation takes 
the original family to a new family g^^ix"" ; e), which satisfies the same requirements as the 
original family. We follow Moncrief [T^ in constructing, from the 3-metric 7^^ on constant-t 
surfaces, quantities g^, which are invariant to first order in e ("gauge invariant"), for coor- 
dinate transformations. The construction of these Moncrief qi is done in two steps. First, 
the multipole moments of the metric are extracted. In practice this is done by multiplying 
the metric functions by certain angular factors and integrating over angles. Since we are 
only interested in quadrupole and higher order for radiation, this step also eliminates the 
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spherically symmetric background parts of the metric function. The second step is to form 
linear combinations of these multipoles and of their derivatives with respect to radius. We 
symbolically represent the process of forming these quantities as 

Qi = Qi{lij,drjij) . (5) 

Here the symbol "Qj" represents the process of multiplying by angular functions and inte- 
grating, then multiplying by certain functions of r and taking linear combinations of the 



results. (Our notation here disagrees with that of Moncrief |jT^ in a potentially confusing 
way. Moncrief 's perturbation quantities are independent of the size of e. In order to have 
definitions that can be applied to numerical data we use quantities that - to first order - 
are proportional to e.) 

The Moncrief gauge invariants play two different roles. For even parity one of the gauge 
invariants, q2, is a constraint; it vanishes in linearized theory as a result of the initial value 
equations. In linearized theory, the remaining Moncrief quantities, denoted qi here, satisfy 
wave equations L{qi) =0, the Regge- Wheeler equation in odd parity and Zerilli equation in 
even parity. 

From our numerical data we construct the quantities qi precisely according to Our 
numerically constructed "perturbation" quantities will not be invariant under coordinate 
transformations, but rather will transform as g- = qi + 0{e'^^j^). Similarly, the linearized 
constraint, q2 will not vanish, but will be of order e^^^. The numerically constructed wave- 
functions qi will satisfy L{qi) = 0((enum)^)j where L is the Regge- Wheeler or Zerilli wave 
operators. 

The use of NFS methods is equivalent to ignoring the second order terms in the wave 
equations. The wavefunction can then be propagated from the initial hypersurface forward 
and the radiation waveforms extracted from it. To evolve qi off the initial hypersurface, 
however, requires the initial time derivative dqi/dt. This can be computed from the initial 
extrinsic curvature, but some care is needed. Indeed, the possible ambiguities that arise 
here are the justification for the somewhat protracted discussion in this section. 

If n is the future-directed unit normal to the initial hypersurface then the rate at which 
the 3-metric is changing is given by 

Kij = -^C^-fij , (6) 

where Kij is the extrinsic curvature and is the Lie derivative along the unit normal. The 
unit normal is related to the derivative with respect to Schwarzschild time by d/dt = S^^"^ n. 
The time derivative of the Moncrief function then can be written 



dq^/dt = S^'^Cn 



= 5i/2£„Qi(7.„97.,/9r) . (7) 

To evaluate the right hand side we need to know how Qi changes when it is Lie dragged 
by n. Since Qi depends only on •jij it might appear that one need only Lie drag •jij to 
find the change in Qi, and that C^Qi = Qi{'Cn'yij,dCnlij/dr). From this it would follow 
that dqi/dt = —2S^^^Qi{Kij, dKij /dr) . It is important to note that this is not the correct 
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relationship between Kij and the cauchy data for the wave equation. The fallacy in this 
procedure lies in the fact that qi must be computed from the 3-metric on a slice for which 
Schwarzschild time is constant (to first order in enum) • Lie dragging by n moves the 3-metric 
to a surface that is not (to first order) a constant time surface. The cure is clearly to 
compare quantities on surfaces of constant t by using Ct = S^^'^Cn- It is the Schwarzschild 
time derivative that commutes with the Schwarzschild radial derivative Ct{d/dr)°' = 0. The 
correct prescription then follows from 

dqi/dt = S'/^C^q, 



= Ql{S'/^Cnl^J,^{S'/^Cnl^,)/^r) 

= -2Q,{S'/'Ki„ d{S'/'K,,)/dr) . (8) 

We note that the perturbed Schwarzschild metric does have a shift vector /5j of order 
e, and in principle the shift vector influences the time development of 'jij according to 
dt'Jij = dt''~fij + 2V(j/3j), where t' is a time coordinate in which the shift vector vanishes. 
But the shift vector can be considered to be "pure gauge." It is necessary if one wants 
a complete specification of the coordinates and the metric components, but its value is a 
matter of choice, and is not necessary for a complete specification of the physics. The initial 
value, and evolution, of the gauge invariant quantity qi is invariant with respect to the choice 
of Pi, and qi carries all the (physically meaningful) information about gravitational waves. 

The evaluation of qi from (^ and dqi/dt from (^ completes the extraction, from the 
numerical data for 'jij, Kij of the cauchy data for the Regge- Wheeler or Zerilli wave equation. 
An alternative procedure arises if one uses the scalar wave-equations derived from the per- 
turbative reduction of the nonlinear wave-equation for the extrinsic curvature which arises 



in a new explicitly hyperbolic form of the Einstein equations [jT^ . In this system, the scalar 
wave equations are one order lower in time derivative from the usual Regge- Wheeler and Zer- 
illi equations, so the Cauchy data consists of the extrinsic curvature and its time-derivative 
(which involves the 3-dimensional Ricci curvature). 

From the above it is clear that linearized evolution should give good accuracy when ap- 
plied to numerically generated initial data with sufficiently small deviations from sphericity. 
For initial data which are known in analytic form one can, of course, apply linearized theory 
even to cases in which initial deviations from sphericity are only marginally small. The re- 
sults in Fig. 2, for example, show that the results of such application of perturbation theory 
give reasonable accuracy for values of fiQ at which an initial horizon is highly distorted. It 
is worrisome to apply linearized evolution to marginally nonspherical initial data, which do 
not, for example, satisfy the constraint ^2 = with reasonable accuracy. Such a procedure 
— linear evolution of nonlinear initial data — has, among other disadvantages, no clear 
theoretical framework. 



C. Calculating radiated energy by "forced linearization" 



We wish to point out here that NFS methods can be used more broadly, and a procedure 
we call "forced linearization" can be applied to numerically generated initial data in a way 
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that amounts to extracting the hnearized part of the data and evolving hnearly. This 
procedure circumvents the difficulty of performing formal linearization to data which is only 
known numerically. We imagine that we start with an initial value problem in which there 
is some adjustable parameter, call it /x, such that /x = corresponds to the Schwarzschild 
initial data. There is no requirement that the family of solutions (7^jy(a;"; /i) be analytic in 
yU at yU ^ 0. There may be additional parameters, call them pj, such as the parameters 
governing the initial momenta of holes. To apply forced linearization we fix the values of the 
Pi and make a choice of /x such that the computed initial data '-)'Jj'^,K™^ are "very nearly 
spherical." One criterion for this would be that q2 is very small. We then interpret this 
initial data as being essentially linearized data, to which the approximation in (^) applies. 
We extract multipoles, form a gauge invariant wave function qi, and evolve it with the Zerilli 
or Regge- Wheeler equation, all as described above. The result of this will be a late-time 
waveform g]''^*^(r, t) and the energy E^^^ that it carries. The next step is to characterize the 
results with a well behaved gauge invariant parameter. To do this we choose some fiducial 
radius r^, and evaluate e^"*^ = gi(rfid, t = 0) the gauge invariant wave function of the initial 
hypersurface at this radius. 

Next, we leave the pi unchanged, but choose a larger value of n for which the numerically 
generated initial data set 74™"^^°', K^^^'^^ is "marginal" in that it corresponds to deviations 
from sphericity large enough so that it differs significantly form linearized initial conditions; 
one sign of this would be that the condition g2 = is significantly violated. For this data set 
we go through the same procedure as above in characterizing the data set by a parameter 
^mrgni ^ ^^(rfid, t = 0). For this marginally spherical initial data we take the solution for the 
wavefunction and energy to be 



The idea underlying this method is that the very nearly spherical data give us the solution for 
for dg^^/de\^=Q. For the marginal initial data set we then need only multiply this initial data 
by the appropriate factor telling us how much larger is the linear part of the nonsphericity 
than that of the very nearly spherical initial data. The success of forced linearization requires 
then that e evaluated at rgd be a well behaved parameterization of the linearized part of 
the nonsphericity in the numerical data. Since our expansion parameter e is the magnitude 
of the perturbation, it will be a good expansion parameter as long as it is evaluated in a 
region where the nonlinear deviations from sphericity are small, i.e., where is a good 
approximation. For this reason it is important that rgd be chosen fairly large. For processes 
of the type pictured in Fig. 1, the deviations from sphericity fall off quickly in radius, so that 
at large enough r one can be certain that the initial data are an excellent approximation to 
linearized data. Evidence for this is that the violations of the ^2 = constraint are always 
confined to small radii. One easily implemented check on the forced linearization procedure 
is to look at the factor gmrs'iyg^n^ and confirm that it is independent of r for r > rgd. In 
Section III we show that this test is easily passed by a numerical example, and that the 
results of forced linearization are essentially the same as those of formal linearized theory. 





8 



III. EXTRACTION OF PERTURBATIONS FROM NUMERICAL DATA 



Here we assume that the reader has numerical solutions for the 3-metric on an approxi- 
mately t=const surface. The first step in applying NPS to numerical results is to transform 
to coordinates which are "almost Schwarzschild" coordinates. It is assumed that the nu- 
merical and Kij are expressed in a coordinate system i?, 9, in which the approximate 
spherical symmetry is manifest. This means that Kij and ratios like 



V76»e ^lee lee 

must be small. They all are, in fact, formally of order enum, so if they are not all reasonably 
small compared to unity there is little reason to think that NPS will work. A Schwarzchild- 
like areal radial coordinate r needs to be introduced. This can be defined as a function of 
Rhy 

r= (^j ieel4>4>dnY^ I^T^. (11) 

where the integral is taken on a surface of constant R. The metric component 7,.^, in terms 
of this quantity, gives us another test of how close the geometry is to that of a constant time 
Schwarzschild slice. The quantity 

r (1 - l/lrr) 

should be nearly equal to the constant 2M, where M is the mass of the spacetime. The 
variability of this quantity in r, 6', and 0, is formally of order enum- 

There are, of course, other ways of specifying the Schwarzschild-like coordinates. We 
could, for example, have defined = All these coordinate choices, however, should 
agree to order enum and are therefore equivalent within a linearized gauge transformation. 

To compute the gauge invariant perturbation functions, we first assume that an Im 
multipole of the 3-metric may be expanded as 

lij = ci(ei),, + 02(62)., + /ii(/i)y + ^{hh + r''KCh)^, + WGCU)^, (12) 

where, for clarity, we have suppressed multipole indices and have replaced Moncrief 's hi and 
/i2 odd parity perturbation functions with ci, C2. The multipole moments Ci,C2,hi,H2,K, 
and G are computed by projection onto the relevant spherical harmonics which can be 
found in Moncrief Explicit formulas for the important special case of even parity, 

axisymmetric perturbations may be found in Ref. ||1U[| . 

For odd parity perturbations, one function can be constructed from the amplitudes ci 
and C2 which is gauge invariant and satisfies the Regge- Wheeler equation (below). 



Q 



(^ + 2)! 

r(^-2)! 



2\dr r 



-• (13) 

r 



The situation for even parity perturbations is more complicated. Two gauge invariant 
functions may be formed out of the multipole moments: 
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ki^K+-{r'^drG-2hi) 
r 



ko — — 
' 2S 



H2 — rdrki 



rS J 



(14) 
(15) 



Prom ki and k2 it is possible to form two new functions, one which is radiative and one 
which is equivalent to the perturbed hamiltonian constraint 

= ArS'^k2 + + l)rki (16) 

g2 = dr[ArS''k2 + 1{£ + l)rk^\ + l{l + l)[2Sk2 + (1 - M/{rS})k^\. (17) 

The scaled function 



Q 



2(£ -!)(£ + 2) 
£(£ + 1) 



(18) 



with 



A = (£- l)(£ + 2) + 6M/r , 

satisfies the Zerilh equation (below). 

The time derivatives of the radiative gauge invariant functions Q^^^ and are found 



by substituting yl — 2M/rKij for in the multipole moment computation and forming 
the same combinations of moments. 

The wavefunctions Q^^ and obey the Regge- Wheeler and Zerilli wave equations 
respectively: 

LQL + ^/QL^O (19) 

^QL + vtQL = (20) 

where the wave operator appropriate to Schwarzschild spacetime is 

92 



(21) 



in terms of the "tortoise coordinate" r* = r + 2Mln(r/2M — 1), and where the potentials 
are given by 

"£(£+1) 6M 



and, 



2M/r) 



= (1 - 2M/r 
1 /72M2 12M 



A2 r5 



(£-l)(£ + 2)(l-3M/r) + 



(22) 

-!){£ + !){£ + 2) 



(23) 

Once the Zerilli and Regge- Wheeler equations are integrated for all the desired £ and m 
modes, the total radiated energy can be calculated from the asymptotic timeseries for Q^^ 
and Qf^: 



dE 
'dt 



■t 00 £ 

E E 



327r 



1=2 m= 



dt 



dt 



(24) 
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IV. EXAMPLE OF PERTURBATION EXTRACTION 



In this section we demonstrate the extraction of a perturbation from a numerical solution 
to the nonlinear constraint equations - the Misner data representing two black holes at a 
moment of time symmetry. The Misner 3-geometry may be written as 

dP = ^(r, e, fio)\S'^dr^ + Pdn^). (25) 

The conformal factor $ is given by 

$(r,^;/io) = l + 2(l + M/2i?)-i J2 ^ ^^(cos^), (26) 

1=2,4,... ^ ^ ^ 

where 

R = (v^+ Vr - 2M)V4 

and 

^ ( 1 V""' V (27) 

' V4E:r=i(smhn/io)-iy' sinhn/io ' ^ ^ 

For this exercise, we pretend that the initial geometry is known only numerically, so 
no explicit formal linearization can be done. The odd parity perturbations vanish in the 
Misner solution. We compute the even parity gauge invariant wavefunction for £ = 2 using 
numerical evaluations of (P^ ) - (PTj). Specifically, we compute K and H2 of ( [T^ ) from 

K = H2= [ dm^Y2o . (28) 



All the other moments in (^) vanish for the conformally Schwarzschild metric of (^3|). 
The function Q20 is evaluated at values of r corresponding to the range r^, = — 20M to 

= 50M. The initial value of Q2Q (along with its time-derivative which is zero for the 
Misner time-symmetric initial data) provides initial values for integration of (pO[). At large 
radius, r = lOOM, the value of dQ2o/dt is used in (|2^) to compute the radiated energy. 

First, in Fig. ^ we show the result of directly computing the gauge invariant function 
Q20 from the nonlinear initial data, integrating the Zerilli equation, and computing the 
radiated energy. For small values of fiQ the agreement with the explicitly linearized data of 



Ref. is excellent. At about fiQ ^ 1.2 the agreement breaks down and the qualitative 
behavior becomes dramatically different. It is interesting to note that the apparent horizon 
encompassing both black holes does not exist for /xq > 1.36, close to the dramatic reversal 
in the energy curve. 

In Fig. ^ the violation of the linearized constraint by the nonlinear data is shown as a 
function of radius. We plot the ratio of the constrained gauge invariant function, qi to the 
radiative function q2 scaled in such a way as to compensate for large violation at r = 2M. 
The value of q2 clearly grows much faster than the radiative variable qi as the separation is 
increased. 
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As discussed in Sec. II, it is possible to obtain the results of formal perturbation theory 
directly from the numerical data without ever making reference to the analytic solution. In 
Fig. ^ we demonstrate the application of the forced linearization procedure to the nonlinear 
Misner data for various values of the fiducial radius rgd. For very small values of fiQ, such 
as fit) = 0.5, the geometry outside the event horizon is everywhere well approximated by 
(I) and forced linearization works even for small values of Vfid/M. When /iq is larger than 
around 1.5, on the other hand, the initial geometry near the horizon contains significant 
nonlinear effects, and large values of rgd/M must be used to get results equivalent to those 
of formal linearized theory. 

As Tfid gets large, the results become indistinguishable from those of formal perturbation 
theory reported in Ref. For rgd = 30M the difference in radiated energy for /io = 3.0 is 
less than 10~^%. This high-accuracy equivalence deserves some explanation. In particular, 
why is forced linearization equivalent to formal linearization with expansion parameter ^2? 
Why is that expansion parameter singled out? The equivalence is a result of two features 
of the way in which the linearizations were done: First, both the formal linearization of 
Ref. |]lT|, and the forced linearization results in Fig. |^, use precisely the same coordinates. 
(The forced linearization results, in fact, are not based on initial values that were generated 
by genuinely numerical means. Rather, the closed form solutions for the Misner metric 
functions were used. The "almost-Schwarzschild" coordinates of the forced linearization, 
were precisely the same as the "almost-Schwarzschild" coordinates in Ref. [|lll). Secondly, 
in the "almost-Schwarzschild" coordinate system, the parameter K2 is, to all perturbation 
orders, the coefficient of the dominant nonsphericity at large radius. Forced linearization 
(in the limit of large rfid) results in a parameterization based on a gauge invariant measure 
of nonsphericity at large radius. It therefore must be proportional to ^2 and produce results 
equivalent to those of the formal linearization of Ref. [^, in which K2 was the expansion 
parameter. 

It should be understood that this does not imply that the parameter K2 is physically 
singled out. A first order change in the "almost-Schwarzschild" coordinates will change the 
coefficient of the dominant large-radius nonsphericity. We might, for example, transform 
from the "almost-Schwarzschild" radial coordinate r of (^) to a new coordinate r' = r[l + 
'^2-^2 (cos ^)]. In this case the coefficient of the leading large r' term in the metric will be 
K2 + 0{k2), and the results of forced linearization with the resulting "numerical" data will 



differ, when perturbations are large, from the results in Ref. ||TT[. The forced linearization 
will have induced an expansion parameter different from K2. 

AMA was supported by National Science Foundation grant PHY 93-18152/ASC 93- 
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FIGURES 



FIG. 1. Spacetime regions for coalescence. The "legs of the trousers" represent the world tubes 

of two compact objects before coalescence; region I cannot be considered to be nearly spherical. 
The objects coalesce in region II which is also highly nonspherical, but lies inside a horizon. Region 
III, above the hypersurface and outside the horizon, can be treated as a nearly spherical spacetime. 

FIG. 2. The effect of a change of expansion parameter. Results are given for the energy 
radiated, as a function of fiQ, during the head-on collision of two black holes (Misner initial data). 
The results of numerical relativity are compared with linearized theory for different choices of 
expansion parameters. 

FIG. 3. Radiated energies from nonlinear Misner data. Energies computed by integration of 
the Zerilli equation are compared for initial perturbations calculated by explicit linearization of 
the Misner data (solid line) and initial perturbations extracted directly from the nonlinear Misner 
data (dashed line). 

FIG. 4. Violation of the linearized hamiltonian constraint. The ratio of gauge invariant func- 
tions qi I q2 scaled by the factor r — 2M is plotted as a function of tortoise coordinate r* . Gurvcs are 
shown for jiq = 0.2, 0.4, 0.6, 0.8, 1.0, 1.2, 1.4. The largest constraint violation occurs for /xq = 1.4. 

FIG. 5. Radiated energies from forced linearization procedure. Radiated energy is plotted as a 
function of Misner separation parameter jiQ for various values of rgd. The curve for rgd = 30M is 
indistinguishable from the formal perturbation theory result of Price and Pullin. 
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